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Abstract 

The van der Waals friction between two semi-infinite solids, a small 
neutral particle and semi-infinite solid is reconsidered on the basis of 
thermal quantum field theory in Matsubara formulation. The friction 
calculated in this approach is in agreement with the friction calculated 
in the framework of dynamical modification of the Lifshitz theory with 
use of the fluctuation-dissipation theorem. This solves the problem 
about the applicability of the Lifshitz theory to the dynamical situa- 
tion. In quantum field theory the calculation of the friction to linear 
order in the sliding velocity is reduced to the finding of the equilib- 
rium Green functions. Thus this approach can be extended for bodies 
with complex geometry. We show that the van der Waals friction can 
be measured in non-contact friction experiment using state-of-the-art 
equipment. 
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1 Introduction 



A great deal of attention has been devoted to the problem of non-contact 
friction between nanostructures, including, for example, the frictional drag 
force between electrons in two-dimensional quantum wells [H El Ej, and 
the friction force between an atomic force microscope tip and a substrate 
mi3000H[ISHIl[H[in]- This is because of it importance for ultra- 
sensitive force detection experiments. The ability to detect small forces is 
inextricably linked to friction via the fluctuation-dissipation theorem. Ac- 
cording to this theorem, the random force that makes a small particle jitter 
would also cause friction if the particle was dragged through the medium. For 
example, the detection of single spins by magnetic resonance force microscopy 
[Tlj , which has been proposed for three-dimensional atomic imaging ^3] and 
quantum computation [16J, will require force fluctuations (and consequently 
the friction) to be reduced to unprecedented levels. In addition, the search 
for quantum gravitation effects at short length scale ^7], and future mea- 
surements of the Casimir and van der Waals forces ^B|, m ay eventually be 
limited by non-contact friction effects. 

In non-contact friction the bodies are separated by a potential barrier 
thick enough to prevent electrons or other particles with a finite rest mass 
from tunneling across it, but allowing interaction via the long-range elec- 
tromagnetic field, which is always present in the gap between bodies. The 
presence of an inhomogeneous tip-sample electric fields is difficult to avoid, 
even under the best experimental conditions 6J. For example, even if both 
the tip and the sample were metallic single crystals, the tip would still have 
corners, and more than one crystallographic plane exposed. The presence of 
atomic steps, adsorbates, and other defects will also contribute to the spa- 
tial variation of the surface potential. This is referred to as "patch effect" . 
The surface potential can also be easily changed by applying a voltage be- 
tween the tip and the sample. An inhomogeneous electric field can also be 
created by charged defects embedded in a dielectric sample. The relative 
motion of the charged bodies will produce friction which will be denoted as 
the electrostatic friction. 

The electromagnetic field can also be created by the fluctuating current 
density due to thermal and quantum fluctuations inside the bodies. This 
fluctuating electromagnetic field gives rise to the well-known long-range at- 
tractive van der Waals interaction between two bodies [T^l 120], and is re- 
sponsible for radiative heat transfer. If the bodies are in relative motion, the 
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same fluctuating electromagnetic field will give rise to a friction which will 
be denoted as the van der Waals friction. 

The origin of the van der Waals friction is closely connected with the van 
der Waals interaction. The van der Waals interaction arises when an atom or 
molecule spontaneously develops an electric dipole moment due to quantum 
fluctuations. The short-lived atomic polarity can induce a dipole moment 
in a neighboring atom or molecule some distance away. The same is true 
for extended media, where thermal and quantum fluctuation of the current 
density in one body induces a current density in other body; the interaction 
between these current densities is the origin of the van der Waals interaction. 
When two bodies are in relative motion, the induced current will lag slightly 
behind the fluctuating current inducing it, and this is the origin of the van 
der Waals friction. 

The van der Waals interaction is mostly determined by exchange by vir- 
tual photons between the bodies (connected with quantum fluctuations), and 
does not vanish even at zero temperature. On the other hand, the van der 
Waals friction, at least to lowest order of perturbation theory, and to linear 
order in the sliding velocity, is determined by exchange of real photons, and 
vanishes at zero temperature. 

To clarify the origin of the van der Waals friction let us consider two 
flat parallel surfaces, separated by a sufficiently wide vacuum gap, which 
prevents electrons from tunneling across it. If the surfaces are in relative 
motion (velocity v) a frictional stress will act between them. This frictional 
stress is related with an asymmetry of the reflection coefficient along the 
direction of motion; see Fig.l. If one body emits radiation, then in the rest 
reference frame of the second body these waves are Doppler shifted which 
will result in different reflection coefficients. The same is true for radiation 
emitted by the second body. The exchange of "Doppler shifted photons" is 
the origin of van der Waals friction. 

From the point of view of the quantum mechanics the van der Waals 
friction originates from two types of processes: (a) Photons are created in 
each body with opposite momentum and the frequencies of these photons 
are connected by vq x = uj\ + u>2, where q x is the momentum transfer, (b) A 
photon is annihilated in one body and created in another. The first process 

(a) is possible even at zero temperature, and it gives rise to a friction force 
which depends cubically on sliding velocity fU\ 122]- The second process 

(b) is possible only at finite temperatures, and gives rise to a friction which 
depends linearly on the sliding velocity. Thus, process (b) will give the main 
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contribution to the friction at sufficiently high temperatures, and at not too 
large velocities. 

In contrast to the van der Waals interaction, for which theory is well 
established, the field of van der Waals friction is still controversial. As an ex- 
ample, different authors have studied the van der Waals friction between two 
flat surfaces in parallel relative motion using different methods, and obtained 
results which are in sharp contradiction to each other. The first calculation of 
van der Waals friction was done by Teodorovich [22] • Teodorovich assumed 
that the force of friction can be calculated as the ordinary van der Waals force 
between bodies at rest, whose dielectric function depend on the velocity due 
to the Doppler shift. However, from the dynamical modification of the Lif- 
shitz's theory it follows [22! that it is not true. Later the same approach 
was used by Mahanty [21] to calculate the friction between molecules. Both 
theories predict wrong nonzero friction at absolute zero of temperature to 
linear order in the sliding velocity. The same nonzero linear friction at zero 
temperature was predicted in [23l2ni- In |2Zj it was shown that the basic 
equation in (2H1 12S] is incorrect, and a correct treatment gives a vanishing 
linear friction at T = OK. Schaich and Harris developed a theory [2E! which 
describes the dependence of friction on the temperature T and on the sep- 
aration d. However in their calculations they made the series of unphysical 
approximations, and for the Drude model their final formula for the friction 
for parallel relative motion gives a divergent result. The friction obtained 
in |3UJ |in] vanishes in the limit of infinite light velocity c — > oo. How- 
ever, at least for short distances, one can neglect by retardation effects when 
calculating the van der Waals friction, as well as van der Waals interaction. 
Pendry J2IJ assumed zero temperature and neglected retardation effects, in 
which case the friction depends cubically on the velocity. Persson and Zhang 
j22] obtained the formula for friction in the limit of small velocities and finite 
temperature using a simple quantum mechanical approach, again neglecting 
retardation effects. In Ref . E7| Volokitin and Persson developed a theory of 
the van der Waals friction based on the dynamical modification of the well 
known Lifshitz theory [19] of van der Waals interaction. In the nonretarded 
limit and for zero temperature this theory agrees with the results of Pendry 
|21j . Similarly, in the nonretarded limit and for small sliding velocity this 
theory agrees with the study of Persson and Zhang [221 • in PlITU] the theory 
was extended to two flat surfaces in normal relative motion. For the case of 
resonant photon tunneling between surface localized states, normal motion 
results in drastic difference from parallel relative motion. It was shown that 
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the friction may increase by many orders of magnitude when the surfaces 
are covered by adsorbates, or can support low-frequency surface plasmons. 
In this case the friction is determined by resonant photon tunneling between 
adsorbate vibrational modes, or surface plasmon modes. When one of the 
bodies is sufficiently rarefied, this theory gives the friction between a flat 
surface and a small particle, which in the nonretarded limit agrees with the 
results of Tomassone and Widom (221 • A theory of the van der Waals friction 
between a small particle and flat surface, which takes into account screening, 
nonlocal optic effects, and retardation effects, was developed in |2I] . 

Dorofeyev et. al. jl] claim that the non-contact friction observed in [3J |3] 
is due to Ohmic losses mediated by the fluctuating electromagnetic field. 
This claim is controversial, however, since the van der Waals friction for 
good conductors like copper has been shown [221 123 E3| to be many or- 
ders of magnitude smaller than the friction observed by Dorofeyev et.al. In 
[35] it was proposed that the non-contact friction observed in Ref.jH] can be 
explained by the van der Waals friction between the high resistivity mica 
substrate and silica tip. However in the experiment the mica substrate and 
silica tip were coated by gold films thick enough to completely screen the 
electro dynamic interaction between the underlying dielectrics. 

At small separation d ~ lnm, resonant photon tunneling between adsor- 
bate vibrational modes on the tip and the sample may increase the van der 
Waals friction by seven orders of magnitude in comparison with the good 
conductors with clean surfaces US]. However, the distance dependence 
(~ l/d 6 ) is stronger than observed experimentally jH]. 

Recently, a theory of non-contact friction was suggested where the friction 
arises from Ohmic losses associated with the electromagnetic field created by 
moving charges induced on the surface of the atomic force microscope tip by 
the bias voltage or spatial variation of the surface potential ^] . It was 
shown that the electrostatic friction can be greatly enhanced if there is an 
incommensurate adsorbate layer which can exhibit acoustic vibrations. This 
theory gives a tentative explanation for the experimental non-contact friction 
data jHj. 

Although at present there are many theories of the van der Waals friction, 
which are frequently contradict to each other, the rigorous theory based 
on the quantum field theory is still not developed. For the van der Waals 
interaction such quantum field theory was developed in Ref.j2ni- On the 
base of this theory were solved some problems which can not be solved in the 
Lifshitz theory of the van der Waals interaction. Due to of great importance 
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of the van der Waals friction for the understanding of the origin of the non- 
contact friction there is strong demand for the creation of the rigorous theory 
of the van der Waals friction. 

This article is organized as follows. In Section 2 we present a short 
overview of the basic idea of the quantum field theory of the van der Waals 
friction. We apply this theory for the calculation of the van der Waals friction 
between two semi-infinite solids (Sec. 3), and a small particle and semi-infinite 
solid (Sec. 4), for parallel and normal relative motion. These calculations 
confirm early results obtained with use of the dynamical modification of the 
Lifshitz theory and the fluctuation-dissipation theorem. Thus these calcu- 
lations solves the problem of applicability of the Lifshitz theory for bodies 
in relative motion. The quantum field theory of the van der Waals friction 
is more general and can be applied for bodies with complex geometry. In 
Sec. 5 we show that the van der Waals friction can be greatly enhanced for 
high-resistivity metals, dielectrics with strong absorption in low-frequency 
region, and for two- dimensional systems, e.g. 2D-electron systems on the 
dielectric substrate, or incommensurate layers of adsorbed ions exhibiting 
acoustic vibrations. The origin of this enhancement is related to the fact 
that screening in 2D-systems is much less effective than for 3D-systems. The 
fluctuating charges will induce "image" charges in the 2D-system. Because 
of the finite response time during relative motion of the bodies this "image" 
charge will lag behind the charge which induce it, and this results in a force 
friction acting on the bodies. The weaker screening in 2D-systems will results 
in larger lag of the "image" charge in 2D-systems as as compare 3D-systems, 
and, as the consequence, in larger friction. Sec. 6 presents the conclusions 
and the outlook. Appendix A contains the details of the derivations of the 
Green functions of the electromagnetic field in two plane surface geometry. 
These Green functions are used in Sees. 3 and 4. 



2 General formalism 

There are two approaches to the theories of the van der Waals interaction 
and the van der Waals friction. In the first approach the fluctuating elec- 
tromagnetic field is considered as a classical field which can be calculated 
from Maxwell's equation with the fluctuating current density as the source 
of the field, and with appropriate boundary conditions. This approach was 
used by Lifshitz in the theory of the van der Waals interaction and by 
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Volokitin and Persson for the van der Waals friction [TO] 122] . The calcu- 
lation of the van der Waals friction is more complicated than of the van 
der Waals force because it requires the determination of the electromagnetic 
field between moving boundaries. The solution can be found by writing the 
boundary conditions on the surface of each body in the rest reference frame 
of this body. The relation between the electromagnetic fields in the different 
reference frames is determined by the Lorenz transformation. The advantage 
of this approach is in that, in principle, it can be used for the calculation of 
friction at the arbitrary relative velocities. However, the calculations become 
very complicated for bodies with complex geometry. At present the solutions 
are known for the van der Waals friction between two parallel plane surfaces 
[TUJ |22j , and between a small particle and plane surface |27j . 

In the second approach the electromagnetic field is treated in the frame of 
the quantum field theory jSH] . This approach was used in Ref . [2U] to obtain 
the van der Waals interaction for an arbitrary inhomogeneous medium all 
parts of which are at rest. 

For two bodies in slow uniform relative motion (velocity v) the force 

acting on either body may be written as F = F — T -v, where the adiabatic 
force Fq is independent of v, and T, the so-called friction tensor, is defined 



Here (...) represents a thermal average of the fluctuating force in the equi- 
librium state at fixed separation d between bodies, and F(£) is the force 
operator in the Heisenberg representation. Eq. Q is a consequence of the 
fluctuation-dissipation theorem [37]. For the interaction between a localized 
and an extended system, Eq.(JTJ) has been derived by several authors (Schaich 
1974 [HE], d'Agliano et al 1975 |Ej], Nourtier 1977 [H]) and is also valid for 
two extended systems. In the context of the van der Waals friction Eq.ljT]) 
was used by Schaich and Harris j2S], but their treatment is incomplete. 

In the case of extended systems the fluctuating force operator can be 
expressed through the operator of the stress tensor 



where the integration is over the surface of the one of the bodies and 



by 




(1) 




(2) 




(3) 
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where Ei and Bi are the electric and magnetic induction field operator, re- 
spectively. The calculation of the force-force correlation function can be done 
using the methods of the quantum field theory EI] • Such calculations 
are described in Sec. 3 for two plane parallel surfaces, and in Sec. 4 for a 
small particle and plane surface, for parallel and normal relative motion. 
The advantage of this approach is that it only involves finding of the Green's 
functions of the electromagnetic field for the equilibrium system with fixed 
boundaries. Thus, this approach can be easily extended to bodies with com- 
plex geometry. However, it is restricted to small relative velocities. 



3 Van der Waals friction between two plane 
surface 

3.1 Parallel relative motion 

Assume that the xy-p\ane coincides with one of the surface. For parallel 
relative motion the friction coefficient T\\ = T xx = T yy . Using the methods 
of quantum field theory j^H] the expression for the friction coefficient for 
parallel relative motion can be written in the form 

r„ = hm l m ^>° + ^, (4) 

where G xx is the retarded Green's function determined by 

G*H = -y o dte** (F x (t)F m {0) - F,(0)F,(0) (5) 

where 

F x = J dS z a xz , (6) 

where the surface integral is taken over the surface of the body at z — 0, 

a xz = (E X E Z + E Z E X + B X B Z + B Z B X ) /8tt (7) 

The function G xx can be obtained by analytic continuation in the upper half 
plane of the temperature Green's function G xx determined on the discrete 
set of point iu n = i2irn//3 by the formula 

G xx {iu n ) =-\£ dre^ T (T t F x (t)F x (0)) , (8) 
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where n is an integer and (3 = H/hsT. T T is the time- ordering operator. The 
function G xx {iuj n ) can be calculated using standard techniques of quantum 
field theory HJJ. Thus, the function G xx {iu n ) can be represented through 
the Green's functions of the electromagnetic field 

D EE (r,r',iuj n ) = Ai(r,r',zuv) = ~ dre^ T (t t ^(t)4(0)) , (9) 

where the retarded Green functions Dij(r,r',uj) obey the equations 

(VjVfc - <5 lfc V 2 ) D kj (r,r',u) - (uj/c) 2 J d 3 x"e ik (r, r", u;)D kj (r", r', a;) 
= (4vRu7c 2 )4,<5(r _ r ') (10) 

(V; V' fc - 5 jfc V' 2 ) Afe(r, r', w) - (^/c) 2 1 d 3 x"e kj (r", r', w) Afc(r, r", u) 
= (47r^ 2 /c 2 )5 ij 5(r - r') (11) 

For the plane surface it is convenient to decompose the electromagnetic field 
into s- and p— polarized plane waves. In this representation with q — q/q 
and n = [z x q], where q is the surface component of the wave vector, the 
Green's tensor is given by 

^■EE f cP(J 

D (r, r') = J —^(fiD EE (z, z', q)n + qD EE (z, z', q)q + zDg E (z, z', q)z 
+ zD zq (z, z', q)q + qD qz (z, z' , q^e****"*') (12) 



where we have taken into account that D^ z = D^ E = (see Appendix A). 
For two plane surface geometry the solution of Eqs. (fTU| [TT]) is described in 
Appendix A. Using the methods of the quantum field theory [SHI E] for the 
Green function G xx we get 

, . TtJ^ f d ^ — ^ Q r , r EE f i\ t~\ EE/ i \ 

G xx {lUJ n ) = -j^p j (2^)2 ^2 [ D qq ^' Z )-°" (~ q ' % ^n~^m, Z, Z ) 

+D EE (q, iu m , z, z)D EE (-q, iu n - iu m , z, z') 



+Df q B (q, iu m , z, z')Df E (-q, iu n - iu m , z, z') 



+ Df z B (q, iu m , z, z')Df q B (-q, iu n - iu m , z, z') 
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(13) 

z=z'=0 



where A is the surface and Df- ( 



(q, ico m , z, z') is given by (HE] 



rot ifc rot^Z>^ E (r, r', zu; n ) 



(14) 



In Eq. (|13|) we omitted terms involving product of Green's functions associ- 
ated with the p- and s- polarized electromagnetic field because after sum- 
mation they cancel to each other. Even without any detailed calculations it 
is clear that such terms must give zero contribution to the friction because 
the p- and s- polarized waves must give independent contributions to the 
friction. 

All the sums over uj m in Eq.()13|) can be calculated in a similar way. Thus, 
as illustration, we consider only one sum: 



According to the Lehmann representation, the Green's function can be writ- 
ten in the form 



pgf (w, r, r') = ex P( F - E n ){E a {v)) nm (E a {v')) mn {l - e -^)8(cu - u mn ) 



Using (|16|) and standard rules for the evaluation of the sum like ()15|) [4*T] we 

get 




(15) 




(16) 



where 



n,m 




/°° r / \ 

du (pf/(q, u)Df E (-d, iu n - u)) n(u 
-oo L v ' 

+ (D EE (q^u n - u)(% E (-q,uj) (n(u) + 1)" 



(17) 



Using Eqs.(jnnj and (|6T|) in (JT7J) we get 



\ E D t E ^ "m)D EE (-q, iu n - iu m ) = 



10 



V- 



du 



7T J-oo 



d p qq E {uo,z,z') d D qq (iuo n -uo,z,z' 



dz> 



+ 



dz< 



7 2 (co>) dz ^ 2 {iuo n — uo) 

n{u) + 1) 



n(uj) 



8 Df q E (lU n -U,Z,Z') d pf q E {uJ, Z, Z 1 ) 



dz 



7 2 M 



;is) 



T{iu n - uo) 

where 7 2 (tu) = {uj/c) 2 — q 2 . Replacing %uo n — > uo + iS and taking the imaginary 
part of Eq. ffTHj) gives, in the limit uo — > 0, the following contribution to the 
friction coming from Eq.()18|) 



limlml. limjj;^^)^ 



2hq 2 

7T7 4 



00 I dn 
duo 







duo I \dz 



\m.D qq {uo 



zq v q,^ 

( d 



lUO r , 



dz , ^D qq {u) 



(19) 



Performing a similar calculation of the other sums in Eq fjlHj) gives 



^7T 3 JO 



00 / <9n" 



o? 2 q g : 



(27r) 2 g 2 



ImD (?(? ImL) 22 



•7 (J^, 



4 2 r # 2 

\mD nn - - - ImD 

nil 



q 



' dzdz' 



'd_ 

dz 1 



lmD r 



~} z =z'=v 



(20) 



Using Eqs. ()66I67I71I73|) for Green's functions in Eq. (j20j) . the contribution 
to the friction from the propagating (q < uo/c) waves becomes: 



rad 



■>7T 3 JO 



Re 



duo I — — 

1 + Ri P R2 P e 2lld — Rip — R 2 pe 2%ld 



1 + Ri p R 2p e 2t,yd + R\p + R 2p e 2tld 



1 - e 2 ^ d R lp R 2p J ' 

Rip — R 2p e 2t,yd 



e 2 ^ d R lp R 2p 



h 

S7T 2 JO 



Im 



00 / dn 



1 - e 2i ^ d R lp R 2p 

ui/c 



+ [p 







dqq c 



(1- 




Rip 


2 )(1- 


R 


2p 


2 ) 




1 


- e 2i ^ d R lp R 2p 


2 



+ [p 



(21) 
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Similarly, the contribution to the friction from the evanescent electromagnetic 
waves (q > uj/c): 

duj I — 



7|1 : 8tt 3 Ju 



{- d £)L,/ ml 



Im ( 2R lp R 2p e~ 2 ^ d - R lp - R 2p e~ 2 ^ d \ ( 2R lp R 2p e- 2 ^ d + R lp + R^e^ 
' m V 1 - e-W d R lp R 2p ) m { 1 - e-^\ d R lp R 2p 

( Ri p - R 2p e~ 2 ^ d \ 2 
ft Z" 00 , / dn\ f°° , o Imi?i„Imi?2iD r i 

= 2^1 ^("a;j/„/ we 171 Rs4r + h!l (22) 

Eqs. fl21l22|) were first derived in Ref. [22] using the dynamical modification of 
the Lifshitz theory with use of the fluctuation-dissipation theorem. 

3.2 Normal relative motion 

For two plane surfaces in normal relative motion the operator of the force is 
given by 

F z = / dS z a zz , (23) 



where 

&zz = {E Z E Z — E X E X — E y Ey + B Z B Z — B X B X — B y B y ) /8n (24) 

The friction coefficient for normal relative motion can be obtained from 
the analytical continuation of the Green function G zz (iu m ) which is deter- 
mined by 



r (it » ^ - ^ A ( V" 9 * [ y\EE T)EE , n EE n EE , r)££ n 



_f) EE f) EE _ f) EE D EE _ f) EB D EB _ f) EB f) EB 

zq zq qz qz zn zn nz nz 

+ Df n B Df n B + D EB D EB + [E^ B)] (25) 

where the arguments of the Green functions in Eq.([25p are the same as in 
Eq. (jl3|) . [E <-> B] denotes the terms which can be obtained from the first 
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terms by permutation of E and B in the upper case indexes of the Green 
functions, and 

(26) 
(27) 



roM il D^{T,r>,u n ) 



rot^ J Dg £ (r,r',^ 



Performing similar calculations as for the parallel relative motion we get 



duo 



dn \ f d 2 q 



x 



16vr 3 y V duJJ ( 2?r ) 2 



' a; 



7 4 (ImD nn ) 2 + 



(. 



2 



\ dzdz' 

2 



lmD r 



+ 2 7 2 (AlmA,„(.,.')) ]}^„ 



(28) 



Substitution the expressions for the Green's functions from Eqs. f)66!67l71l73|) 
in Eq. (|2*K|) gives the contribution to the friction from propagating waves: 



16vr 3 



00 , / dn 



J Jq<ix}/c 



1 + Ri p R2 P e 2l ' yd — R\ v — R2 V e 2lld 



Re- 



+2 Im 



1 - e 2 ^ d R lp R 2p 



Rip — R2r > e 2l ' yd 



1 - e 2 '^ d R^R 



lpJTL2p / 



+ [P 



+ ^Re 

n 



1 + RipR2 P e 2t " ld + + R2 V e 2nd 



l p t -n-2pt 



47T 2 Jo 



v l — \R\ p \ 2 \R2p 



2\2 



+ |(1 - |i? lp | 2 )i? 2p e^ + (1 - \R 2v \ 2 )RLe-^ d \ 2 



1 - e 2 ^ d R lp R 2p 

dn\ t w l c , 2 
dqq'j 



+ [P 



^ U [ Qu) Jo 



2 P \ J-n-ipt 



|l-e 2 ^i? lpJ R 2p | 4 

(29) 

In similar way one can obtain the contribution to the friction from the evanes- 
cent electromagnetic waves: 



_ evan 

7± 



471-3 
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X 



' lm lR lp R 2p e- 2 ^ d - R lp - R 2p e~ 2 ^ d J | ^2R lp R 2p e~ 2 ^ d + R lp + R 2p e~ 2 W 



1 - e-W d Rx p R 2p J V 1 - e- 2 \-<\ d R lp R 2p 

_ 2 (i^'-^'f*) 2 ^ - .]i = \ r<t, 

^ l- e -^\ d R lp R 2p ) 1 JJ Wo V duJJui/c yy|/l 
x[(IiilR 1p + e- 2|7|d| | J R lp | 2 Im J R 2p )(Im J R 2p + e" 2|7|d| |i? 2p | 2 Imi? lp ) 

+ e-H-^^OT—^^ + ^ s] (30) 

Eqs. f|29l30|) were first presented without derivation in Ref.p. The derivation 
was given in Ref. ^U] on the base of the dynamical modification of the semi- 
classical Lifshitz theory ^5] of the van der Waals interaction and the Rytov 
theory [121 HHl HI] of the fluctuating electromagnetic field. 

4 Van der Waals friction between a small par- 
ticle and plane surface obtained using quan- 
tum field theory 

4.1 Parallel relative motion 

Let us assume that the xy- plane coincides with surface of the body and that 
the z- axes is directed along the upward normal. For parallel relative motion 
the friction coefficient F\\ = T xx = T yy . The Lorentz force acting on a small 
particle located at point ro can be written in the form 



k^-E x {x) + - { Jy B z - Jz B y ) (31) 

^ J r=ro 

where p and j are the dipole moment and current operators of the particle, 
respectively. E and B are the external electric and magnetic induction field 
operators, respectively. The interaction of the electromagnetic field with the 
particle is described by the Hamiltonian 

H int = —A(r )-j (32) 
c 

where A(r) is the vector potential operator. Taking into account that 

j = |p (33) 
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V xE=~B 

c at 



(34) 



one can prove that the friction coefficient is determined by Eg. (HI) where 



h Jo 



dte 



iujt 



d d do 

x (p k (t)—E k (v, t^M—E^', 0) -nW—Etf, 0)p k (t)—E h (v, t))^^ 

(35) 

where summation over repeated indexes is assumed. After similar calcula- 
tions as in Sec. 3 we get 



2h 

71 



\ / k=x,y,z 

^1 } r =r'=r 



2Re (a xx (uj)a* zz (uj)) I — ImD K (r, r , u 



(36) 



where -Djj(r, r') are the Green's functions of the electromagnetic field for one 
plane surface geometry without interaction with a particle. These Green 
functions can be obtained from the Green functions for two plane surface 
geometry (see Appendix A) by putting Rp( s ) — 0. The polarizability of the 
particle 

a kk (u) = T / dte tujt < Pk (t)p k (0) - p k (0)p k (t) > (37) 
n Jo 

is expressed through its value a kk in the absence of interaction between the 
particle and the surface 



a«(w) 



1 - a%(uj)Du(r ,ro) 



(38) 



We have also used the identity 
lma xx (uj)lm 



d d 
a zz (u)—D xz (r, r , u)—D xz (r, r , u) 



+Ima 22 (u;)Im 



d d 
a xx (u)—D xz (r, T , u)—D xz (r, r , lo) 



-21m 







acxx(v)-z-D xz (r,r ,u}) 
ox 



Im 







9 

a zz (u)—D xz (r,r ,u) 
ox 



2Re (a xx (u)a zz (u)*) I — lmD xz (r, r , u) 



(39) 
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4.2 Normal relative motion 



The friction coefficient for a particle moving normal to the sample surface 
can be obtained from calculations very similar to those for parallel relative 
motion. In this case the Green's function G% x must be replaced by Gf z where 



dte 



d d d d 

x (p h (t)-^ E k {r, t) Pl (0)—Ei(x', 0) - p z (0)— ^(r', 0) Pk (t)—E k (r, t) 

After similar calculations as in Sec. 4.1 we get 



r=r'=ro 

(40) 



On 



IT JO 



d 2 



k=x,y,z 



dzdz' 



lmD kk (r, r', u) 







+lm(a k kD k k(r,r ,uj)D kk (r',r ,uj)) + —lm(a kk (uj)D kk (v,r ,uj)) \ 

J L(jz J J r=r'=r 

(41) 

For a spherical particle with radius R Eq. (|41j) is only valid if R <C d. In non- 
resonant case a° kk ~ R 3 and D kk ~ d~ 3 . Thus in this case we can neglect by 
the screening effects. Than taking into account that for a spherical particle 
&kk — and using in the non-retarded limit (which can be formally obtained 
as a limit c — > oo) formula (see [22] and also Appendix A) 



D kk (r,r',u) = 47r 

k=x,y,z 



we get 



r,, = 2- / 

7T JO 



d 2 qq 

WT 2 



dn(uj) 
du 



e- qlz - z ' l +R p (q,uj)e- 



-q{z+z') 



3 iq(x-x') 



(42) 



dqq A e 2qd lmR p (q, w)Imal 



u 



(43) 



and = 237(1 . Eqs. (|3l^ particle9) were first obtained in Ref. [2ZJ with use of 
the dynamical modification of the semiclassical Rytov theory ^21 G31 IE| of 
the fluctuating electromagnetic field. The particular case of these equations 
given by Eq. (j4^ was first derived in Ref.[33] using the fluctuation-dissipation 
theorem. 



16 



5 Van der Waals friction between dielectrics 
and two-dimensional systems 

In Refs. [3*4*1 ITU] it was shown that the van der Waals friction between good 
conductors (hsT/Ana ^> 1, where a is the conductivity) is extremely small. 
However the van der Waals friction can be greatly enhanced for high resis- 
tivity materials. Recently [!""""] the calculations were published, where the 
authors claim that they can explain the experimental data about long-range 
noncontact friction observed by Stipe et.al [""]. In this experiment the sub- 
strate did not consist of a bulk conductor but of a metal film with a thickness 
of a few hundred nanometers deposited on top of a dielectric substrate. Thus 
in the theory from Ref. |3Hj it was proposed that experimental measurements 
performed on metal films do not reflect the properties of the metal but of 
the underlying dielectric substrate. Using the formula ("4*3*]) with the particle 
polar azability 

a(u) = R s£ -^ (44) 
v ; e+2 v ; 

and the reflection coefficient in the electrostatic limit (q ^> ui/c) 

Hp~, («) 
for high- resistivity material (Ana << k%T jK) we get: 

9 k B T R 3 1 



2 Ana d 5 2e' + 3 



(46) 



where e' is the real part of the dielectric function e = e' + 4nia/uj. In the 
frequency range below 0.75 x 10 12 Hz the real part of the dielectric constant 
of glass is nearly constant [e' ~ 3.82), and a = 30s -1 . It was assumed 
that d — s + h s + h t , where s is the tip's apex-surface separation, and h s 
and ht are the thicknesses of the gold films which were deposited on the 
mica substrate and on the silicon cantilever, respectively. In experiment jH] 
h s = 250nm, h t = 200nm, R = 1/xm, than for s = lOnm and T = 300K we get 
T|| = 7 x lCT 13 kg/s, which agrees exactly with value obtained numerically 
in [S3], and in rough agreement with the experimental value reported in 
Ref. jOj. However the macroscopic theory which was used in obtaining Eq. ()46j) 
is only valid when the average separation between electrons is much smaller 
than length scale of variation of the electric field, which is determined by 
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separation d. Thus the lowest value of the electron concentration n min is 
restricted by the condition n min > d~ 3 and according to Drude formula the 
lowest value of the conductivity a > a m i n ~ e 2 r/md 3 . For d ~ 1/xm and 
t = 10~ 15 the conductivity must be at least four orders of magnitude larger 
than it was used in [33] ■ For such conductivity Eq.pfij) gives T» ~ 10 _16 kg/s, 
which is three orders of the magnitude smaller than the observed friction. 

Besides, in the experiment jB] the mica substrate and silica tip were coated 
by gold films thick enough to completely screen the electrodynamic inter- 
action between underlying dielectrics. If planar film with thickness h and 
dielectric function 63(00) lies on a substrate with dielectric function £2(^0, 
the system response is similar to the single-interface case if one replace the 
reflection coefficients R P 2i by 

R = R P 3i - i? P 32exp(-2g/i) 

p l-R p31 R p32 exp(-2qh) 1 ] 

where 

R pij = ~ ~ ~ ( 48 ) 

Sj ~\~ Si 

. The magnitude of the wave vector q is determined by the separation between 
metallic films. Thus qh ~ h/s ^> 1, and in this case reflection coefficient 
Rp « R 3l instead of R21 as it was proposed in Ref . |35j . 

For high-resistivity metals (hsT/Ankcr > 1), for d < c(h/4ira/lll the 
reflection coefficient (??) in Eq. (|22jl for two surfaces in parallel relative we 
get 

7 ~ 0.05^^ (49) 

and 7j_ ~ 107||. From the condition of the validity of the macroscopic theory 
(see above) maximum of friction can be estimated as 

h k B T mk B T 
7||ma * ~ °- ^4^- ~ 4^d- (50) 

To estimate the friction coefficient T for an atomic force microscope tip 
with the radius of curvature R » d we can use the "proximity approxima- 
tion" |HH 147] . Thus the maximal friction coefficient for a spherical tip: 

q ma ,~o.i 7mfl ,^~^^. (si) 
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For r ~ l(r 16 s, R ~ 1/im and T = 300K we get T max ~ l(T 15 kg/s. This 
friction is two order of magnitude smaller than the friction observed in a 
recent experiment at d = lOnm [6 . Similarly, in the case of a cylindrical tip: 

-rc fTB mk B TR l/2 w 

T maa; ~ Ima^dRw ~ ^2^1/2 ( 52 ) 

where w is the width of the tip. For w = 7/xm, d = lOnm, and with the 
other parameters as above, the friction is of the same orders of magnitude as 
it was observed in the experiment 0. 

Recently a large electrostatic non-contact friction between an atomic force 
microscope tip and thin dielectric films was observed jT3|. Similarly the van 
der Waals friction will be also large between dielectrics with high absorption 
in low-frequency range. As a particular important case we consider the van 
der Waals friction between thin water films adsorbed on the surfaces of the 
transparent dielectric substrates like silica or mica. Water has an extremely 
large static dielectric function of around 80. The low frequency contribution 
to the dielectric function, responsible for this large static value, is due to 
relaxation of the permanent dipoles of the water molecules. It is very nicely 
described by the simple Debye [48J rotation relaxation. The theoretical fit of 
the experimental data is given by 



e'{u) = 4.35 + — - — (53) 

(1 + {uJ/UJqY) 

C{u/u ) 2 

£ M = (54) 

where C = 72.24 and ljq = 1.3 • 10 11 s _1 . The fit is very good for both the real 
and imaginary part of the dielectric function e{oj) = e'{uj) + ie"{uS) We note 
that water has large absorption in the radio-frequency range at u ~ ujq, and 
shows in this region of the spectrum anomalous dispersion. In this frequency 
range the dielectric constants e<i of mica or silica are nearly uniform and 
I £3 1 3> £2, where £3 denotes the dielectric function of water. The reflection 
coefficient for substrate with film is given by Eq. (j47|) . For qh 1 and 
q^ 1 ~ d <^ \£$\h/e2 it can be approximated by equation 

R p ^l (55) 

e 3 qh 
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Substituting in eq.(j22J) and using "proximity" approximation for friction 
between cylindrical atomic force microscope tip and a sample we get 



irhR l ' 2 w (k B T 



Q^2C 2 h 2 d?/ 2 V hu 



(56) 



For h = lnm and with the same parameters as above we get rjj = 4.8 • 
10 _12 kg/s what is one order of magnitude larger than the friction observed 
in Ref. [U]. It is interesting to note that the friction coefficient (J56|) has the 
same weak distance dependence as it was observed in Ref.jB]. 

Another enhancement mechanism of the van der Waals friction is con- 
nected with resonant photon tunneling between adsorbate vibrational modes 
localized on different surfaces. In [HI GUI we have shown that resonant pho- 
ton tunneling between two surfaces separated by d — lnm, and covered by 
a low concentration of potassium atoms, result in a friction which is six or- 
ders of the magnitude larger than for clean surfaces. The adsorbate induced 
enhancement of the van der Waals friction is even larger for Cs adsorption 
on Cu(100). In this case, even at low coverage {6 ~ 0.1), the adsorbed layer 
exhibits an acoustic branch for vibrations parallel to the surface (35]. Thus, 
uj\\ = and according to Ref. ^2] at small frequencies the reflection coefficient 
is given by 

R p = l- - \ q (57) 
or — u A q + za>?7 

where uj 2 = 4im a e* 2 aq 2 /M, e* is the ion charge and a is the separation 
between an ion and the image plane. Substituting Eq. (|57jl in Eq. (|22|) for 



a A7m a e* 2 a 

Md 2 ^ 

and using a "proximity" approximation, for a cylindrical tip we get 

T^0. e S kBTa T W (58) 

II ^5.5 V > 

For Cs adsorption on Cu(100) the damping parameter 7] was estimated in 
[TJ] as r] w 10 u s _1 . Using this value of rj in Eq.(|55jl for a = 2.94A |33], 
R = lpm, w = 7/im, T = 293 K at d = lOnm we get T {1 w -lO" 15 kg/s, 
which is two orders of magnitude smaller than the friction observed in jH] 
at the same distance. However, the van der Waals friction is characterized 
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by a much stronger distance dependence (~ 1/d 5 ' 5 ) than observed in the 
experiment (~ l/c? n , where n = 1.3 ± 0.2). Thus, at small distances the van 
der Waals friction will be much larger than friction observed in JBJ, and can 
thus be measured experimentally. 

6 Summary and conclusion 

We have used a quantum field theory in Matsubara formulation to calculate 
the van der Waals friction between two plane surfaces, and a small parti- 
cle and plane surface, for parallel and normal relative motion. The friction 
calculated in this approach is in agreement with the friction calculated in 
the framework of dynamical modification of the Lifshitz theory with use of 
the fluctuation-dissipation theorem. This solves the problem of the appli- 
cability of the Lifshitz theory to the dynamical situation. In quantum field 
theory the calculation of the friction to linear order in the sliding velocity 
is reduced to the finding of the equilibrium Green functions which obey to 
the system of the Maxwell type differential equations. Thus with applica- 
tion of the numerical methods of classical electrodynamics this approach can 
be used in the calculations of the van der Waals friction between bodies 
with complex geometry. We show that the van der Waals friction between 
high-resistivity metals, dielectrics with strong absorption in radio-frequency 
range, and two-dimensional systems can be measured in non-contact fric- 
tion experiment using state-of-the-art equipment. The theory can be used as 
a guide for planning and interpretation of new experiments including more 
sophisticated dielectrics like in Ref.|T3]. 
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A The Green functions for two plane surface 
geometry 

Suppose that the half-space z < is occupied by a solid with the reflec- 
tion coefficients -Ri p (q, uj) and Ri s (q,cu) for s- and p- electromagnetic fields, 
respectively. Similarly, the half-space z > d is occupied with a solid with 
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the reflection coefficients R 2p (q,,uj) and R 2s <l,u)). The region < z < d is 
assumed to be vacuum. Here q is the surfaces component of wave vector 
k = (q, 7), where 7 = ((to/c) 2 — q 2 ) 1 ^ 2 . Since the system is uniform in the 
x = (x, y) directions, the Green function Dy(r, r') can be represented by the 
Fourier integral 



1 2 

A;(r,r', me<7<) = J —^e^^D 



yZ, z , q, uv n ) 



(59) 



In the xy— plane it is convenient to choose the coordinate axes along the 
vectors q = q/g and n = z x q. In this coordinate system the equations (jlOl 
ITT|l for the Green functions become 



7 2 + 



dz 2 



D (z z' 



UJ 



dz 2 



a 



D qq (z,z') - iq—D zq (z,z' 
d 

-y 2 D zq (z, z') - iq—D qq (z, z 
d 

-i 2 D zz (z, z') - iq—D qz (z, z' 
7 2 D qz (z, z) + iq—D qq (z, z 



Attlo 



c 

Attu 2 



2 5 i z ~ z ') 
2 5{z-z') 



5(z - z') 



(60) 
(61) 
(62) 
(63) 
(64) 



The components D qn , D zn of the Green function vanish, since the equations 
for them turn out to be homogeneous. Solving the system of Eqs. (j60|) - (}64j) 
amounts to solving two equations: equation (|HUJ) for D nn and the equation 
for D qq which follows from equations (|HT| E2J 



7 2 + 



8T 
dz 2 



D qq (z,z') = -4n 1 2 5(z-z>), 



(65) 



after which D qz , D zq and D zz for z 7^ z' are obtained as 



D R 

qz 



■ d = i±JL D . 

-fdz' qv zq 7 2 <9z «' 
q 2 d 2 
-fdzdz' qq 



(66) 
(67) 



22 



w n = R 2s e 2 ^ d ( v n + —^e-^ for z = d (70) 



In the vacuum gap < z < d the solution of equation ()6Ujl has the form 

D nn {z, z') = Jj^-^-A + Vn6 ^ + Wrifr *i* ( 6 8) 

At the boundaries z = and 2 = d the amplitude of the scattered wave is 
equal to amplitude of incident wave times to corresponding reflection coeffi- 
cient. The Green function D nn is associated with s- polarized electromagnetic 
field and the boundary conditions for it give 

v n = Rn(w n + ^^^A for z = (69) 

/ 2-rriuj 2 x 

■jc 

Using Eas. (IH%H?0l we get 

D nn (z,z>) = J^L^-A + R ^ e I 6 +e J 

7c 2 { A s 

+ A s J (71) 

A, = l-e 2 ^ d R 2s R ls (72) 

Equation for D qq is similar to equation (JoTIj) for -D nn , and the expression 
for D^ q can be obtained from expression (J71|) just by replacements of reflection 
coefficient 

Dw = © Dnn [Rs ~* ~ Rp] (73) 

In our approach the calculation of reflection coefficient for s- and p- polar- 
ized waves constitutes separate problem, which can be solved taking into 
account non-local effects. For the local optic case the reflection coefficients 
are determined by the well known Fresnel formulas 

Rip ! j Ris ! j 

where Si is the complex dielectric constant for body i: 



7i = v-^ £ i-i 2 ( 75 ) 
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FIGURE CAPTION 

Fig. 1. The electromagnetic waves emitted in the opposite direction by 
the body at the bottom will experience opposite Doppler shift in the reference 
frame in which the body at the top is at rest. Due to the frequency dispersion 
of the reflection coefficient these electromagnetic waves will reflect differently 
from the surface of the body at the top, which give rises to momentum 
transfer between the bodies. This momentum transfer is the origin of the 
van der Waals friction. 
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